Draft version March 13, 2008 

Preprint typeset using I^T^X style cmulateapj v. 10/09/06 



oo 
O 
O 
(N 



43 

6 



> 

OS 



o 

oo 
o 



X 



THE SPECTRAL LAG OF GRB 060505: A LIKELY MEMBER OF THE LONG DURATION CLASS 

S. McBreen 1 , S. Foley 2 , D. Watson 3 , L. Hanlon 2 , D. Malesani 3 , J. P. U. Fynbo 3 , D. A. Kann 4 , N. Gehrels 5 , 

S. McGlynn 2 , D. Palmer 6 
Draft version March 13, 2008 

ABSTRACT 

Two long 7-ray bursts, GRB 060505 and GRB 060614, occurred in nearby galaxies at redshifts of 
0.089 and 0.125 respectively. Due to their proximity and durations, deep follow-up campaigns to 
search for supernovae (SNe) were initiated. However none were found in either case, to limits more 
than two orders of magnitude fainter than the prototypical GRB-associated SN, 1998bw. It was 
suggested that the bursts, in spite of their durations (~ 4 and 102 s), belonged to the population of 
short GRBs which has been shown to be unrelated to SNe. In the case of GRB 060614 this argument 
was based on a number of indicators, including the negligible spectral lag, which is consistent with 
that of short bursts. GRB 060505 has a shorter duration, but no spectral lag was measured. We 
present the spectral lag measurements of GRB 060505 using Suzaku's Wide Area Monitor and the 
Swift Burst Alert Telescope. We find that the lag is 0.36 ± 0.05 s, inconsistent with the lags of short 
bursts and consistent with the properties of long bursts and SN-GRBs. These results support the 
association of GRB 060505 with other low-luminosity GRBs also found in star-forming galaxies and 
indicates that at least some massive stars may die without bright SNe. 
Subject headings: gamma rays: bursts 



1. INTRODUCTION 

The existence of two classes of gamma-ray bursts 
(GRBs) differing in observed durations and spectral 
properties has been established for s o me time (e.g . 
iMazets et all fl98lt iNorris et all Il98l iHurlevI 11991 . 
These populations were quantified using the Burst and 
Transient Source Experiment (BATSE), which showed a 
bimodal distribution in t he durations of GRBs well fit by 
two lognormal fu nctions (McBre en et al.l[1994h . with the 
divide at ~ 2 s (jKouveliotou et al .11 1993T K In addition, 
there is also contamination in the short burst class f rom 
soft gamma-ray repearters (e.g. lChapman et ai1l2008h . It 
is generally accepted that long GRBs have their origins 
in massive star progenitors because of their association 
with core-col l apse supernovae (SNe, [ Galama e t al. I [19911 
Hiorth et all 120031: IStanek et all 120031: [MalesanGLaD 



GRBs would be associated with SNe (|Zeh et al.l [2004 
iWooslev fc Blooml l200ll) . as predicted by the collap - 
sar model of long GRBs ([Wooslev fc MacFadvenlll999h . 
Doubts were cast on this paradigm by the non-detection 
of SN e in two nearby GRBs, GRB 060505 at z = 
0.089 (lOfeket al.1 120061 and GRB 060614 at z = 0.125 
(iPrice et al.ll2006T ) discovered by the S wift (IGehrels et al.l 
[2004J) Burst Alert Telescope (BAT, iBarthelmvl [20051 . 
Due to their lo ng durations, Tgn of 4 ± 1 s and 102 ± 
5 s re spectively (jHullinger et al.l 120061 : iBarthelmv et al.1 
2006), SN searches were initiated. Although a supernova 
~100 times fainter than SN1998bw w ould have been de- 
tected, none was found in either case (iFvnbo et al1l2006l: 
Gal- Yam et~aHl2006t iDella Valle eFa l. 2006; Ofc k et all 



20071) . It was suggested that they were short bursts 



2004t lPian et al. 2006; Wo oslev k Bloom 2006) and oc 



currence in star- forming galaxies ( Bloom et al.ll2002l) and 
in hi ghly star-forming regions therein ([Fruchter et al.l 
2006). The origin of short GRBs is still open, with merg- 
ers of compact obi ects being the leading concept (e. g . 
IGehrels et all [20051 [Hiorth et al.l[2005l : iFox et alJl2"005h . 

The detection of the spectroscopic signatures of SNe 
in the 4 nearest GRBs and the detection of bumps con- 
sistent with SNe in the lightcurves of most low-redshift 
bursts seemed to confirm the paradigm that all long 
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where the lack of SNe would not be surprising, as short 
GRBs have not shown SN emission (Hior th et all 1 2005: 
IFox et al.ll2005t iBloom et al.ll200l ICovino et al.ll2006fL 

The classification of GRBs with durations close to the 
long/short division is problematic. The argument that 
GRB 060614 was a "short GRB" rests on its extended 
soft emission c ompo nent and on i ts neg ligible spectral lag 
([Gehrels et al.ll2006l : IZhang et al.ll2007f l When the latter 
is combined with its relatively low l uminosity, it violates 
the lag-luminosity relation found bv lNorris et all ((2000a) 
for long GRBs. If the lack of a SN in GRB 060505 is to be 
attributed to it being a short burst, it should also have 
a negligible lag. We present the spectral lag analysis of 
the prompt emission of GRB 060505 using data from the 
Suzaku Wide Area Monitor (WAM) and Swift-BAT. 

2. OBSERVATIONS AND DATA REDUCTION 

GRB 060505 was detected by the BAT instrument on 
Swift. The fluence is (6.2±1.1) x 10 -7 ergs cm -2 (15- 
150 keV) and the spectrum is fit b y a power law with 
index 1.3±0.3 ([Hullinger et al.ll2006[ ). The trigger fell be- 
low the 6.5cr threshold for an automatic slew but ground 
analysis found a 8.5tr excess (|Palmer et al.l[2006h . Swift 
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Fig. 1.— The lightcurves of GRB 060505 with the BAT instru- 
ment on Swift (a) and from the WAMO detector on Suzaku (b-e), 
all at 100 ms resolution. Time is relative to the BAT trigger. A 
precursor is visible in the lowest energy WAM channel. 

was repointed at Tp+ 0.6 days and a weak fad ing X-ray 
source was identified (jConciatore et al.ll2006fh We ob- 
tained the publicly available data for GRB 060505 from 
the Swift archive 7 . A mask weighted lightcurve was gen- 
erated using the BAT data analysis tools. The available 
data contained only 10 s of event data and the lightcurve 
is presented in Fig. [TJ Swift was approaching the South 
Atlantic Anomaly when the burst occurred and was sub- 
ject to a higher than normal background level. Addition- 
ally, the partial coding was only 11% (Hullingc r et al.l 
2006) meaning that the off-axis angle with respect to the 
Swift axis was almost 50° , substantially reducing the ef- 
fective area of the instrument. Splitting the data into 
energy channels for spectral lag analysis further reduces 
the weak signal. 

The WAM is the anti-coincide nce shield (ACS) of the 
Hard X-ray Detector on Suzaku (jYamaoka et al.l l2006a) 
and it also triggered on GRB 060505. The WAM consists 
of four identical walls which act as individual detectors 
(WA MO to WAM3). The de tectors have a large effective 
area (jYamaoka et al.ll2006ah . They are sensitive in the 
energy range 50-5000 keV, and although its primary role 
is to act as an ACS, WAM is also used as an all-sky mon- 
itor for G RBs. An automated tri ggering system operates 
on board ( Yamaoka et al. 2006b) and the lightcurves are 
publicly available at 15.6 ms resolution in 4 rough en- 



ergy bands 50-110, 110-240, 240-520, 520-5000 keV 8 . 
The lightcurves in the four energy channels from the 
WAMO detector at 100 ms resolution are presented in 
Fig. [TJ The T 90 of GRB 060505 was -4.8 s in the 50- 
5000 keV band 8 . The burst struck the WAM detector at 
an angle such that principally WAMO, but also to a lesser 
extent WAM3, detected the burst. The on-axis effective 
area of t he BAT and WAM ins truments are shown in 
Fig. 2 of lYamaoka et al.l (|2006af ) and the effective area 
of WAM only exceeds that of BAT above 300 keV. How- 
ever, it should be remembered that GRB 060505 occurred 
~ 50° off-axis in BAT and that the effective area of BAT 
also drops rapidly above 100 keV. These factors result in 
a more significant detection GRB 060505 by WAM than 
BAT and therefore we rely primarily on the WAM data 
for our analysis. However, we show that the results are 
consistent with those obtained from the BAT data. 

3. DATA ANALYSIS AND RESULTS 

The spectral lag was calculated by cross-corre latin 
the l i ghtcurves in different energy channel s (e.g. iBan 
119971 iNorris et al.ll2000al : iFolev et al.ll2007[ L The cross 
correlation function (CCF) was fit with a fourth order 
polynomial and the quoted lag value is the peak of this 
function. Statistical errors w ere calculated using a boot- 
strap method as described in lNorris et al.l |2000a). This 
involves adding Poissonian noise based on the observed 
counts to the lightcurves in the different energy channels 
and re-computing the CCF in 100 realisations for each 
burst. The 50th ranked value is the mean lag and the 
16th and 84th ranked values represent ±ler. 

The spectral lag was determined between the 50-110 
and 110-240 keV (rno-240, 50-110) energy bands for the 
Suzaku WAM detectors over a range of temporal res- 
olutions (15.6, 31.2, 46.8, 62.4, 78 and 100ms). The 
lightcurves were correlated from —4 to +4s and the CCF 
was fit over a range of —5 s. A lightcurve threshold 
of 10% (30%) is applied, which means that only data 
with at least one-tenth (three-tenths) of the peak count 
rate is used to calculate the lag, thus reducing the back- 
ground. The spectral lag values obtained from WAMO at 
the six time resolutions specified above at 10% threshold 
agree within the uncertainties and the average value is 
0.36 ± 0.05 s. Applying the 30% threshold to the same 
lightcurves, increases the average value to 0.44±0.06s 
and all values are within ler of those obtained at the 10% 
threshold except at 60 ms and 78 ms, which are consis- 
tent at the 2cr level. Above the 50% threshold the results 
are unreliable and the lag is not accurately reproduced. 
The burst is detected with lower significance in WAM3 
and does not allow an accurate determination of the lag. 
We add the signal from WAM3 to that of WAMO to test 
if this gives a consistent result. The results are consistent 
with WAMO alone within — la, except at 100 ms resolu- 
tion where the WAMO+3 lag is larger but is consistent 
at - 3cr with the WAMO results (10% threshold). The 
average value obtained from the sum of the WAMO and 
WAM3 lightcurves is 0.42±0.05s and 0A7t°m s at 10% 
and 30% thresholds respectively. The cross correlation 
data and fit for WAMO at 100 ms is presented in Fig. [2]a 
and is inconsistent with the negligible lag expected for 
a short burst. A precursor is evident in the WAM data 



7 http:/ /swift. gsfc.nasa.gov/docs/swift/archive/ 
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Fig. 2. — The cross-correlation data and fit for a) the WAMO data 
at 100 ms between the 110-240 and 50-110 keV energy bands and b) 
the BAT data at 100 ms between the 50-100 and 25-50 keV bands. 
A 4 th order polynomial fit to the data is shown. The vertical lines 
denote zero lag. GRB 060505 is clearly inconsistent with zero lag. 

at — 8s and including it in the lag analysis over a wider 
time range results in a consistent lag measurement of 
0.47±0.06s. We note that precursors are not normally 
detected in short bursts. 

The lag was also measured between the 25-50 keV and 
50-100keV energy bands (t5 _i 00j 25-50) at 100ms us- 
ing the BAT data with the techniques outlined above 
(Fig. H]b). The lightcurve was too weak to determine 
the lag at finer time resolution. The spectral lag value 
of 0.4±0.1s measured using the BAT data is consis- 
tent within 1 a with that obtained from the WAMO and 
WAMO+3. 

In order to establish the robustness of our result, we de- 
termined the lag for 16 additional GRBs detected by both 
BAT and WAM, for which the lightcurve data were suf- 
ficient for lag analysis in both instruments. The analysis 
was performed in a similar manner to GRB 060505. The 
derived lags ranged from — 3 ms to 0.94 s in the WAM and 
to 0.86 s in the BAT. The lags are compatible consider- 
ing the differing instruments and off-axis angles, energy 
ranges and the spectra of the bursts, except in 2 cases 
where the BAT lag was significantly longer. The sample 
consisted of 12 long and 4 short bursts. Crucially, the 
short bursts were always found to have negligible lag in 
both instruments. This shows that our analysis is sensi- 
tive to short lags. 

4. DISCUSSION 

4.1. Spectral lags 

The spectral lags in G RBs have b e en discussed 
by many authors (e.g. I Band! Il997t iNorris et al.l 
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Fig. 3 . — The lag-luminos ity relation using data (diamonds) and 
fit from INorris et al.1 l!2000aD . In addition, GRB 060 505 (open cir- 
cle), GRB 060614 (filled square: iGehrels et al.ll2006l ). short GRBs 
(open squares) and 3 GRBs associated with SNe (filled-circles) are 
included. The la g values are fr om the following: GRB 060218: 
[ Liang et al]||2006h . GRB 031203: IFolev et all p)07h . GRB 980425: 
INorris! [120021 ). The spectral lag of GRB 060505 is significantly 
longer than those measured for short GRBs, and it falls on the lag- 
luminosity plot in a position similar to that of some SN-GRBs. The 
diamond and filled-circle lag values are determined between the 
25-50 and 100-300 keV energy ranges. Lags for GRB 060614 and 
the short bursts are measured between the 15-25 and 50-100 keV 
ranges. No fc-correction is applied. 



20008^: INorris k BonnelH 1200(1 .lHakkila k Giblinl [20061: 
Hakkila et alj|2007t IFolev et al.ll2007l) 

Using BATS E data, INorris! (|2002f) and 

INorris &; Bonnell (|2006l ) found that long duration 
GRBs had both measurable and ze ro lags but that 
short GRBs had lags around zero. INorris &: Bonnell 
(2006) calculated the lags of 260 short GRBs using 
BATSE data and found that 90-95% of the values 
were consistent with zero and suggest that bursts with 
positive lag may result from contamination by the long 
GRB class. It was also argued that if short GRBs had 
lags proportionally as large as long GRBs, such lags 
would be detectable, i.e. that this result was not simply 
an effect of the duration of short bursts. This is not to 
say that bursts with short lags are necessarily in the 
short GRB class. In the sam ple of published lags of 
BATSE data by lHakkila et all (|2007f) . 1427 bursts have 
Tqq > 2s and a measured lag (tioo-5o,20-5o)- Of these 
bursts 214 have lags in the range from —10 to +10 ms 
(79 with uncertainties of ±10 ms) and 348 have lags in 
the range from —20 to +20 ms (217 with uncertainties of 
±20 ms), showing that there are many long GRBs with 
very short lag. In summary, long bursts are expected 
to have predominantly positive lags ranging from zero 
to several seconds. Short GRBs have almost exclusively 
negligible lags. However, it is not possible to exclude 
that GRB 060505 could be an outlier: i.e. a short 
duration GRB with a positive lag or due to a process 
which does not fit into the lag classification scheme. 

There have been difficulties in classifying a number 
of bursts and the lag has been used to discri minate 
in a number of cases (e.g. iDonaghv et alj 12006). For 
example, GRB060912A has a T90 of ~ 6 s and was 
initially thought to have occurr e d in a nearby ellipti- 
cal galaxy, however iLevan et al.1 (|2007l ) recently found 
that it was more likely to come from a star forming 
galaxy at z = 0.937 and report a lag (tioo-350, 25-50) 
of 83 ± 43 ms. Various strategies have been proposed 
to distinguish bursts more effectively than the dura- 
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tion alone (eg. INorris et al.l l2000at IDonaghv et"alll2006t 
IZhang et al.ll2007t ). However, none have seen widespread 
adoption. 

4.2. What was the progenitor of GRB 060505? 

It was argued that GRB 060505 was probably part of 
the tail of the short burst population and connected to 
mergers of compact objects. At a redshift of z=0.089, 
GRB 060505 has an isotropic peak luminosity of ~ 9 x 
10 47 ergs _1 (50-300 keV). Having a low luminosity and 
relatively long lag of 0.36±0.05s, GRB 060505 f alls be- 
low the lag-luminosity relation of INorris et alj ([2000a) 
as shown in Fig. [3l The spectral lag of GRB 060505 is 
significantly longer than those measured for short GRBs 
and GRB 060614 and it falls on the lag-luminosity plot in 
a position similar to that of some (but not all) SN-GRBs 
( e.g. GRB 03 1 203). 

lOfek et all (120071) argue that the simplest interpreta- 
tion for GRB 060505 is that it is related to a merger event 
rather than a short-lived massive star and point out that 
the maximum allowable distance of GRB 060505 from a 
star- formin g knot is cons i stent with the shortest merger 
timescales. rThone et all (|2008l ) claim that GRB060505 
occurred in a star forming region of the host galaxy which 
resembles long GRB host galaxies and argue for a mas- 
sive star origin for this event. It has also been argued 
that the host galaxy of GRB 060505 is more similar a 
short burst host in terms of m etallicity and ionisation 
state (|Levesque fc Kewlevl 12007). However, their short 
GRB host region in the emission line ratio diagram is 
based on only two burst host galaxies, one of which is 
the GRB 050416 A, which has photome tric evidence for 
an associated SN (|Soderberg et aLlfeOOTl ) and is argued to 
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be a long GRB due its spectral softne ss and location on 
the Amati plot (jSakamoto et al.ll2006h . The host galaxy 
studies alone do not resolve the classification issue for 
GRB 060505. The optical luminosity at 12 hours in the 
source frame is similar to those of short GRB afterglows, 
but optical luminosity alone i s also not a valid classifi- 
cation tool (|Kann etal] 12001 . In our opinion, the lag 
measurement suggests that this burst is similar to long 
GRBs implying a massive star progenitor, despite the 
lack of a SN detection. 

It has been argued that the absence of a SN signa- 
ture in GRB 060505 is evidence of a new, quiet endpoint 
for so me massive stars (jFynbo et al.ll2006t IWatson et al.1 
I2007D . The existence of a SN was a feature of the early 
collapsar model. However, the complete absence of a 
SN may be expected where the 56 Ni does not have suf- 
ficien t impetus to escap e the black hole ()Heger et al.1 
120031 : I Fryer et ail l2006f) or in jet-induced explosions 
with narrow jets when the deposited ener gy is small 
(jNomoto et al.l l2007t iTominaga et~aT1 120071 ). Progeni- 
tor stars with relatively low angular mom entum could 
also produce GRBs without supernovae (|MacFadvenl 
2003). These seem attractive explanations at least 
for GRB 060505. In the absence of a GRB explosion, 
the detection of such quiet-death massive stars, if they 
exist, is a challenge for current instrumentation (e.g. 
Kochanek et~alll2008l ). 
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